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Abstract 



We study the curvaton scenario in supersymmetric framework paying particu- 
lar attention to the fact that scalar fields are inevitably complex in supersymmetric 
^ ' theories. If there are more than one scalar fields associated with the curvaton mecha- 

. nism, isocurvature (entropy) fluctuations between those fields in general arise, which 

may significantly affect the properties of the cosmic density fluctuations. We exam- 
ine several candidates for the curvaton in the supersymmetric framework, such as 
moduli fields, Affleck-Dine field, F- and D-Hat directions, and right-handed sneu- 
trino. We estimate how the isocurvature fiuctuations generated in each case affect 
the cosmic microwave background angular power spectrum. With the use of the 
recent observational result of the WMAP, stringent constraints on the models are 
derived and, in particular, it is seen that large fraction of the parameter space is 
excluded if the Affleck-Dine field plays the role of the curvaton field. Natural and 
well-motivated candidates of the curvaton are also listed. 



1 Introduction 



Study of the origin of the cosmic density fluctuations is a very important subject in cos- 
mology. In the conventional scenarios, inflation is assumed as a mechanism to provide 
the source of the cosmic density fluctuations. In the inflationary scenarios, the universe is 
assumed to experience the epoch of de Sitter expansion in the early stage. During the de 
Sitter expansion, physical scale expands faster than the horizon scale so the homogeneity 
of the universe is realized at the classical level. At the quantum level, however, scalar 
field responsible for the inflation, called "inflaton," acquires quantum fluctuation, which 
becomes the origin of the cosmic density fluctuations. If inflation provides the source of 
the cosmic density fluctuations, energy scale of the inflation is related to the amplitude of 
the cosmic density fluctuations. Importantly, observations of the cosmic microwave back- 
ground (CMB) anisotropy give us important informations about the inflationary models. 
In particular, in the simplest scenarios, inflation models should reproduce the currently 
observed size of the CMB anisotropy of AT/T ~ O(10~^). Furthermore, spectral index for 
the scalar-mode perturbations should be close to 1 so that the shape of the CMB angular 
power spectrum becomes consistent with observations. After the very precise measure- 
ment of the CMB angular power spectrum by the Wilkinson Microwave Anisotropy Probe 
(WMAP) jl], thus we obtain stringent constraints on the inflation models. 

Recently, a new mechanism is proposed where a late-decaying scalar condensation 
becomes the dominant source of the cosmic density fluctuations 013111 El- In this scenario, 
a scalar field, called "curvaton," other than the inflaton, acquires primordial fluctuations 
in the early universe. Although the energy density of the curvaton is subdominant in the 
early stage, it eventually dominates the universe. Consequently, the isocurvature (entropy) 
fluctuations originally stored in the curvaton sector become the adiabatic ones ^ which 
become the dominant source of the density fluctuations of the universe. Then, when the 
curvaton decays, the universe is reheated and the fluctuation of the curvaton produces the 
adiabatic density fluctuations. 

The curvaton mechanism has important implications to particle cosmology. If the 
curvaton mechanism is implemented in inflation models, the size and the scale dependence 
of the cosmo logical density fluctuations become different from the conventional result. In 
particular, the curvaton mechanism provides a natural scenario of generating (almost) 
scale-invariant cosmic density fluctuations which is strongly suggested by observations. 
As a result, we can relax the observational constraints on the inflation models 

Various particle-physics candidates of the curvaton field have been discussed so far. 
(For the recent discussions on the curvaton scenario, see Ref. 0.) Among them, it is often 
the case that the curvaton mechanism is considered in the framework of supersymmetry, 
since the supersymmetry can protect the flatness of the curvaton potential which is re- 
quired for a successful curvaton scenario. A crucial point in supersymmetric theories is 
that scalar fields are inevitably complex and hence they have two (independent) degrees 
of freedom. Thus, if the curvaton mechanism is implemented in supersymmetric models, 
effects of two fields should be carefully taken into account. In particular, if there are two 
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scalar fields, isocurvature fluctuations between them are in general induced which may 
signiflcantly affect the behavior of the cosmic density fluctuations. 

Thus, in this paper, we study implications of the curvaton scenarios in the framework 
of supersymmetric models. We pay a special attention to the fact that, in the supersym- 
metric framework, curvaton mechanism requires at least two (real) scalar flelds. The CMB 
anisotropy in such a framework is studied in detail and we investigate how the CMB power 
spectrum behaves. 

The organization of the rest of this paper is as follows. In Section|21 we flrst present the 
framework of the curvaton mechanism. Several basic issues concerning the cosmic density 
fluctuations are discussed in Section IHl Then, cases where the cosmological moduli flelds, 
Affleck-Dine fleld, F- and Z^-flat directions, and the right-handed sneutrino play the role 
of the curvaton are discussed in Sections IH IHl and [71 respectively. Section |S1 is devoted 
for the conclusions and discussion. 

2 Framework 

We flrst present the framework. In the scenario we consider there are two classes of 
scalar flelds which play important roles. The flrst one is the inflaton fleld which causes 
inflation.*^ In the conventional scenarios, the inflaton fleld is assumed to be responsible 
for the source of the cosmic density fluctuations but here this is not the case. The second 
one is the curvaton which becomes the origin of the cosmic density fluctuations. Here, we 
consider general cases where there exist more than one curvaton flelds. Hereafter curvaton 
is denoted as (with relevant subscripts as mentioned below). 

In this paper, we consider the case where the universe starts with the inflationary 
epoch. During inflation, the universe exponentially expands and the horizon and flatness 
problems are solved. In addition, the curvaton flelds are assumed to have non-vanishing 
amplitude during inflation. Assuming that the masses of the curvaton flelds are much 
smaller than the expansion rate of the universe during inflation, the energy density of the 
curvatons are minor component in the very early universe. 

After inflation, the universe is reheated by the decay of the inflaton fleld and radiation 
dominated universe is realized. (We call this epoch as "RDl" epoch.) At the time of 
the reheating, the curvaton flelds are minor components and their energy densities are 
negligibly small compared to that of the radiation. As the universe expands, the expansion 
rate decreases and the curvatons start to oscillate when the Hubble parameter becomes 
comparable to the masses of them. Then, the (averaged) amplitudes decrease. 

In discussing the evolution of the curvatons, it is convenient to decompose the com- 
plex fleld into real ones. One convenient convention is to use the mass-eigenstate basis. 
When the amplitudes of the scalar flelds become small enough, the potential of can be 

Alternative models such as the pre-big-bang [51 [5j and the ekpyrotic ^U] scenarios were proposed. 
However, they predict unwanted spectrum for the density fluctuations Thus, we assume inflation 

in the following discussion. 
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approximated by the parabolic one as F ~ m^|0p + m'^(0^ + h.c). (Notice that m'^ can 
be always chosen to be real using the phase rotation of the field.) Then, the real and 
imaginary parts of 4> become the mass eigenstates. Expanding as 

^=-^(01+^02), (2.1) 

the scalar potential is described by the quadratic form around its minimum as 

V-l^ml^l (2-2) 

i 

where (pi denotes i-th mass eigenstate with i — 1,2. Here and hereafter, the "hat" is 
for complex scalar fields while the scalars without the "hat" are understood to be real. 
In addition, the indices i, j, ■ ■ ■ are for specifying the mass eigenstates. In general, this 
basis should be distinguished from the ones expanding the fiuctuation during the infiation. 
Hereafter, we assume that the masses of the (pi fields are of the same order of magnitude. 

Although the mass-eigenstate basis is useful in particular in discussing the behavior of 
the curvaton fields when the amplitudes become small, it is sometimes more convenient 
to use the "polar-coordinate basis" where is decomposed as 

= ^ (0init + <Pr) e^(^'""+'^«/'^'"'*). (2.3) 

Here, (pmit/V^ and ^init correspond to the initial amplitude and phase of the complex scalar 
field (p, respectively. In addition, (pr and (pQ are real scalar fields. 

Adopting the quadratic potential, the (averaged) pressure of each curvaton vanishes 
once the curvaton field starts to oscillate: 

P<t>. = {l<Pl-lml<pf)^ = 0, (2-4) 

where the "dot" denotes the derivative with respect to time t. This means that, at this 
epoch, the energy density of the curvatons behave as that of non-relativistic matter, namely 
P0. oc with a being the scale factor, while the energy density of the radiation drops 
faster, oc a~^. Thus the curvaton fields can dominate the universe, which is one of 
the indispensable conditions for the curvaton scenario to work. (We call this epoch as 
"0 dominated" or "(/)D" epoch.) Then, the curvaton fields decay and the universe is 
reheated again. Consequently, the universe is dominated by the radiations generated from 
the decay products of the curvatons. (We call the second radiation dominated epoch as 
"RD2" epoch.) 



3 Density Fluctuations and CMB Power Spectrum 

During the inflationary epoch, all the scalar-fleld amplitudes acquire quantum fluctuations. 
If one calculates the two point correlation function of a scalar fleld ip with mass m^: 

(0|M^,^)M^,y)|0)inf = / f\5^{t,k)\'e''^^^-^^, (3.1) 
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where x and y are comoving coordinate, then the Fourier amphtude Sip{t,k) is given by 



where is the expansion rate of the universe during the inflation and k is the comoving 
momentum. For the fluctuations we are interested in, the wavelengths are much longer 
than the horizon scale during inflation. 

As one can see from Eq. (j3.2|) . for the superhorizon mode {i.e., for k <^ aHi^f), fluctu- 
ations of the scalar-field amplitude are suppressed when the scalar mass is comparable to 
or larger than H[^f. In inflation models in the framework of the supergravity, it is often 
the case that the effective mass of O(ifinf) is generated to flat directions. (Such a scalar 
mass is called the "Hubble-induced scalar mass.") If this is the case, the curvaton mech- 
anism does not work. However, even in the supergravity, there are several cases where 
curvaton fields (and other flat directions) do not acquire Hubble-induced masses. One 
case is the so-called D-teim inflation ^2] where the "vacuum energy" during the inflation 
is provided by a D-term potential. In this case, flatness of the curvaton potential is not 
severely disturbed since the Hubble-induced scalar masses originate from the F-term in- 
teraction.*^ Another possibility of maintaining the flatness of the curvaton potential is 
the no-scale type inflation. If the Kahler potential is in the no-scale form, it is known 
that the Hubble-induced mass is suppressed [13]. The no-scale type Kahler potential can 
be naturally realized in the sequestered sector scenario jHj where the inflaton sector lives 
in the hidden brane which is geometrically separated from the observable brane where 
the standard- model particles live.*^ In the conventional inflation scenario, in fact, such 
a scenario is problematic since the decay of the inflaton field reheats the hidden brane 
not the observable brane. In the curvaton scenario, however, this is not a problem since 
dominant part of the matter in the universe [i.e., radiation, cold dark matter (CDM), 
baryon asymmetry, and so on) is generated from the decay of the curvaton fields.*^ 

The curvatons also acquire quantum fluctuations during the inflation. As we mentioned 
before, the curvaton is expected to be a complex field in the supersymmetric framework and 
hence there are two (or more) independent real scalar fields associated with the curvaton 
mechanism. If the flatness of the curvaton potential is not disturbed by the inflation 

After inflation, the inflaton starts osciflation, which induces F-term potential. Thus, the Hubble- 
induced mass term appears even for D-term inflation, and the curvaton fluctuations decrease as a~^/^ 
until reheating. Importantly, the amplitude itself is also proportional to a~^/^ with the parabolic potential. 
Thus, the ratio S(p/(j) becomes a constant of time and the effects of the Hubble-induced mass term do not 
change the following discussion. 

'^^■^Such a set-up may be naturally flt into the scenario of the anomaly-mediated supersymmetry breaking 
|14[ 1151 since the hidden and observable branes are usually introduced in the anomaly- mediated 
supersymmetry breaking models to explain the smallness of the Kahler-induced supersymmetry breaking 
scalar masses. 

scenarios with late-time entropy production, in fact, this is a general property and the inflaton 
field does not have to decay into the standard-model particles. Thus, in this case, the inflaton can be a 
hidden-sector field and its interaction with the standard- model particles can be absent. 




(3.2) 
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dynamics, all the curvatons are expected to acquire the quantum fluctuations as given in 
Eq. ()3.2p with <^ -f^inf- In the following, to make our discussion clearer, we consider 
the case where there are two independent real scalar fields, corresponding to the real and 
imaginary parts of the complex scalar field. (Our formulae can be easily generalized to 
the cases with more scalar fields.) 

During inflation, we expand the curvaton (p around its amplitude [(p]inf- 

4> = Minf + 50. (3.3) 

In general, there are two (uncorrelated) real scalar fields parameterizing S(j). We denote the 
fluctuations of such scalar fields as ScpA- (Hereafter, the indices A, B, ■ ■ ■ are for specifying 
the basis expanding the scalar field during inflation.) 

One way of parameterizing 50 is to use the mass basis. When the amplitude of the 
scalar field is small enough, the potential can be well approximated by the parabolic one 
as Eq. ()2.2j) . which enables us to use this basis in order to parameterize 50. One should 
note that, if the scalar potential has interaction terms, primordial fluctuations in 0j may 
generate fluctuations in 0^ with i ^ j. Of course, if the potential is completely parabolic, 
this does not happen. 

Importantly, fluctuations of two independent scalar fields are uncorrelated. Thus, in 
calculating the two point correlation function of the cosmic temperature fluctuations from 
which the CMB anisotropy is obtained, we can regard (Fourier modes of) the fluctuations 
of the curvatons as independent statistical variables, leading to 



(50^505) 



2n 



Sab. (3.4) 



These fluctuations become the origin of the cosmic density fluctuations. If 6(f) a and 50^ are 
independent, the resultant CMB angular power spectrum is the same for any choice of the 
basis. As we will see below, they may provide both adiabatic and isocurvature fluctuations 
and hence a large class of scenarios are severely constrained from the observations of the 
CMB angular power spectrum. 

Once the fluctuations of the curvaton amplitudes are generated as in Eq. (j3.4|) . the 
density and metric perturbations are induced from [50] inf. To parameterize the density 
fluctuations, we define the variable*^ 

Sx ^ ^, (3.5) 
Px 

where the subscript X is to distinguish various components. Hereafter, X = 7, c, b, and 
m are for the photon (or more precisely, relativistic matter), CDM, baryon, and total non- 
relativistic matter, respectively. The density fluctuation 6px is defined in the Newtonian 



"^^We adopt the notation and convention used in Ref. ^7] unless otherwise mentioned. 
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gauge where the perturbed hne element is given using the metric perturbations \1/ and $ 

as 



-(1 + 2^)dr2 + (1 + 2^)6ijdx'dx^ 



(3.6) 



Evolutions of the metric and density perturbations are governed by the Einstein and 
Boltzmann equations. The Einstein equation provides the relation between metric and 
density perturbations: 



1 



2M? 



a^PT 



St + + ujt)Vt 

k 



(3.7) 



and the equation relating $, \E', and the anisotropic stress perturbation of the total matter 



A;2(^ + $) 



1 



(3. 



where ~ 2.4 x 10^*^ GeV is the reduced Planck scale. Here, the subscript "T" is for the 
total matter and the variables px and Vx are the pressure and the velocity perturbation of 
the component X, respectively. In addition, ut = pt/pt is the equation-of-state parameter 
for the total matter, and 



n 



a 
a 



(3.9) 



where the "prime" denotes the derivative with respect to the conformal time r. 

In order to discuss the CMB angular power spectrum associated with the fluctuations 
of the curvaton fields, we first evaluate the density fluctuations at the deep RD2 epoch. In 
such an epoch, the mean free path of the radiation is very short and hence the fluctuation 
of the radiation component becomes locally isotropic. In this case, the anisotropic stress 
perturbation vanishes and we obtain 



V 

7 



4 

— fcK, 
3 ^ 



4$', 



-k6r + km. 



(3.10) 
(3.11) 



In addition, if the anisotropic stress perturbation vanishes, $ = — ^ due to Eq. ()3.8|) . We 
use this relation to eliminate For non-relativistic component, the Boltzmann equation 
becomes 



v 

m 



-kVm - 3$', 

-nv^ + km. 



(3.12) 
(3.13) 
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It is convenient to use the mass basis in discussing the density fluctuations in the RD2 
epoch. Furthermore, it is useful to define the primordial entropy between i-th curvaton 
and the radiation generated from the decay of the inflaton field: 



Si 



(3.14) 



where S^.^^ is the density fluctuation of the radiation generated from the inflaton fleld. 

Notice that, in studying effects of the primordial fluctuations of the curvaton fields, 5p-y.^j 
vanishes in the RDl epoch and hence Si is determined by the initial value of the density 
fluctuation of the curvatons. Since the entropy between and 7inf is conserved for the 
superhorizon mode. Si becomes the entropy between the components generated from the 
decay product of and those not from (p. Then, we deflne the transition matrix T as 



Si = [50a] inf- 



(3.15) 



The transition matrix T depends on models and initial conditions. If the potentials of the 

scalar fields are well approximated by the parabolic one, it is convenient to choose mass 
eigenstate as S4>a- In this case, with relevant unitary transformation, T can be expressed 



as 



= [diag(20r\ 202 \ • • OIm, 



or equivalently 



Si 



4>i 



(3.16) 



(3.17) 



inf 



In this case, the correlations between the fluctuations of different mass eigenstates vanish. 
For general scalar potential, however, T is not guaranteed to be diagonal in this basis and 
the correlation function {6(f)i6(f)j) may become non- vanishing for i ^ j. Consequently, the 
isocurvature ffuctuations for different mass eigenstates can become correlated. 

Now we study the evolution of the density and metric fluctuations with non-vanishing 
Si with Sj — ioi j i. (Note that the effects of individual S^ can be treated separately as 
long as the linear perturbation theory is valid.) For this purpose, it is convenient to expand 
the fluctuations as functions of kr since we are interested in behaviors of superhorizon 
modes. In addition, it is important to note that, since the curvaton fields are minor 
components in the RDl epoch, the total density fiuctuation at that epoch is negligibly 
small and hence [^^'^'^'*^]rdi vanishes. (Here and hereafter, the superscript (50^) is for 
quantities generated from the primordial fluctuation of 0^^.) The density fluctuation of 
the photon generated from the decay product of 0j (j ^ i) is obtained by using the fact 
that the velocities of the scalar-fleld condensations are higher order in kr relative to ^ 
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and 5^.*^ Thus, from Eqs. (imH) and (pnT|) . we obtain 5^^f^^ = 4^(^<^^) + 0{k'^T'^) for 
i 7^ where 7^ denotes the photon generated from the decay products of 0^- Also, the 
density fluctuation of the photon generated from (pi is obtained using the relations [iZj 

[5t]rd2 = -2[^]rd2, (3.18) 

and 

[5t]rD2 = /7J'^7JrD2, (3.19) 
k 

where f^^ = [P'y,. / P'y]RD2 is the energy fraction of the radiation generated from the decay 
products of k-th mass eigenstate. If Si is the only source of the cosmic density fluctuations, 
S, = 1(5^^ - S^^) for J ^ i, and hence we obtain = -f^^f'^^^ 

In general, the primordial fluctuation in 0^ generates fluctuations in various mass eigen- 
states and hence, in the RD2 epoch, the metric perturbation is related to the isocurvature 
fluctuations as 

[^(''^-)]RD2 = -^E/7,^f^"\ (3.20) 

where 'if^^'l'^^ is the metric perturbation generated from [50A]inf. In particular, when i = 1 
and 2 with /^^ + = 1; "^-e can write 



[^]rD2 = -I 



(3.21) 



where 



5 = ^ (5i + ^2) , 5i2 = S^- S2. (3.22) 

Notice that 5*12 becomes the entropy between components generated from the decay prod- 
ucts of (pi and 02. 

In discussing the CMB anisotropy, the CMB angular power spectrum is usually used: 
{AT{x,^)AT{x,^')),= ^Y.i'^l + l)QPi{j-j'), (3.23) 

with AT(x, 7) being the temperature fluctuation of the CMB radiation pointing to the 
direction 7 at the position x and Pi is the Legendre polynomial. In the curvaton scenario. 



T^^Note that it is possible to define those photons generated from the decay products of each curvaton 
and infiaton separately, as far as we are concerned with the density fluctuation in the long wavelength 
limit. The reason for this is that the density fluctuation over the horizon scale is determined only by 
the gravitational potential, as can be seen from Eq. (|3.10() . If we take into consideration the effect of 
the higher order in fcr, such a distinction becomes meaningless since those photons constitute a single 
component fluid with a common velocity perturbation. 
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the cosmic density fluctuations originate from tlie primordial fluctuations of the curvaton 
fields. In particular, as mentioned before, more than one curvaton fields are expected in 
the supersymmetric case. Then, denoting the contribution of dcpA as Cf'^^\ the CMB 
angular power spectrum has the form 

Ci = ^[^('''^■^f [C/^'^-")]|^|=i. (3.24) 

A 

Notice that Ci^'^'^^ oc 6(f)\. Since the fluctuations of different curvatons are uncorrelated 
as shown in Eq. ()3.4|) . no term proportional to 6(j)A^4'B with A ^ B is expected. 

The CMB angular power spectrum depends on the isocurvature fluctuations in the 
baryon and the CDM as well as the metric perturbation. Although the isocurvature 
fluctuations can be separately defined for the CDM and for the baryonic component, the 
shape of C/^"^^^ is, up to normalization, determined by the ratio 

r c{5(t>A) 

'^m — 

where Smr is given by 



(3.25) 

RD2 



Smr — — /T'^'T" (3.26) 



Notice that the isocurvature fluctuations are in general correlated with the metric per- 
turbations. Thus, the form of the CMB angular power spectrum for one curvaton case is 
written in the form 



^1 -j^ / 771 L TTl I \ / 



_ y(adi) _|_ 9^(<59iA)r'('=°rr) _,_ ^(5</'A)2(^(iso) 

Here, Ci^'^^^ and Ci^^°^ agree with the angular power spectra for purely adiabatic and isocur- 
vature fluctuations, respectively, while C^'^"^^^ parameterizes the effects of the correlation. 
Their formal definitions are given as 

1, |S|=0 ' 
0, |S|=1 ' 

_ ^(adi) _ Miso)\ (o r,Q\ 

V/S=l, |*| = |5|=1 ) ■ l^O.ZOJ 

When two (or more) curvatons exist, the above equation is generalized as 

Ci = (^[m^'^-)f^ + 2«:,o,,C^'^'^) + t^soCt^) ■ (3.29) 

Combining this equation with the Eq. (j3.24p . /tcorr and /tiso are related to k^^*^-*) as 



^(adi) _ 


[Cl]\q,\ = 


^Oso) _ 




^(corr) _ 





(3.30) 
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Figure 1: The CMB angular power spectra with correlated mixture of the adiabatic and 
isocurvature fluctuations. Here, we take Km — oo, 10, 3, 1, 0, — 1, — 3, — 10, — oo. 
The overall normalizations are taken as [l{l + l)Ci/27r];=io = 1. 
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Figure 2: Probability contours for mixture of adiabatic and correlated isocurvature fluc- 
tuations. Note that we consider only the region of kiso > |/«corr|- 

Note that Kiso > |«^corr| is always satisfied from their definitions. 

In Fig. ^ we plot the CMB power spectrum Ci for several values of Km with the relation 
/tcorr = i^iso = i^m- (In our study, we use the cosmological parameters flbh'^ = 0.024, 
^mh"^ = 0.14, h = 0.72, and r = 0.166, where Qb and flm are density parameters of 
baryon and non-relativistic matter, respectively, h the Hubble constant in units of 100 
km/sec/Mpc, and r the reionization optical depth, which are suggested from the WMAP 
experiment |18,. We also neglect the scale dependence of Si^'^^\ which is expected to be 
small in a large class of inflationary models.) As one can see, for non- vanishing value 
of Km, the shape of the CMB angular power spectrum may significantly deviate from the 
result with the adiabatic fluctuation. Since the observed CMB angular power spectrum by 
the WMAP experiment is highly consistent with the prediction from the adiabatic density 
fluctuation, a stringent constraint on kiso and Kcorr can be obtained. 

As one can expect from Fig. Q if k^so or /tcorr becomes too large, the CMB angular 
power spectrum deviates from the adiabatic result. In order to derive constraints on these 
parameters, we perform likelihood analysis using the WMAP data. In our analysis, we 
first calculate the CMB angular power spectrum for given values of kiso and Kcorr- Then 
we calculate the likelihood variable using the numerical program provided by the WMAP 
collaboration with the WMAP data Assuming that the probability distribution is 
proportional to the likelihood variable, we obtain the probability distribution on the Kcorr 
vs. Kiso plane, where the total probability integrated over the plane with the uniform 
measure is normalized to be unity (Bayesian analysis |2()j). 
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The probability contours are shown in Fig. |21 From the figure, it is seen that Hiso = 0.4 
is allowed at 95 % C.L. for the pure (uncorrelated) isocurvature case. We can also see that 
a larger value of Kiso is allowed if positive correlation (/tcorr > 0) exists. In the following 
analysis, we derive constraints on various curvaton scenarios using the result presented in 



4 Cosmological Moduli Fields as Curvaton 

We first consider the case where the moduli fields, which are fiat directions in the super- 
string theory, play the role of the curvaton. Interactions of the moduli fields are expected 
to be proportional to inverse power (s) of the gravitational scale, M*, and their poten- 
tial is lifted only by the effect of the supersymmetry breaking. If the moduli fields have 
non- vanishing amplitude in the early universe, they might dominate the universe at a 
later stage and reheat the universe when they decay. Originally, such moduli fields are 
considered to be cosmologically dangerous since their lifetimes may be so long that they 
might decay after the big-bang nucleosynthesis (BBN). In this case, the decay products 
of the moduli fields would spoil the success of the BBN . However, it was pointed out 
that such a problem may be avoided if the lifetimes of the moduli fields are shorter than 
~ 1 sec This is the case if the masses of the moduli fields are heavier than about 
10 TeV or so. If such cosmological moduli fields exist, they play the role of the curvaton 
if the primordial fiuctuations are generated in the amplitudes of the moduli fields during 
infiation. 

Let us first consider the isocurvature fiuctuation in the baryonic component. In dis- 
cussing the effects of the primordial fluctuation in the curvaton, we here neglect the pri- 
mordial fluctuation in the baryonic sector, which may become independent and uncorre- 
lated source of the cosmic density fluctuation. In order to parameterize the isocurvature 
fluctuation in the baryonic sector, we deflne 



Fluctuation in the baryonic sector strongly depends on when the baryon asymmetry 
of the universe is generated. If the baryon asymmetry is generated after the epoch is 
realized, primordial fluctuation in the curvaton is imprinted into the baryonic sector and 
the isocurvature fluctuation in the baryonic component does not arise. This happens when 
the baryon asymmetry is somehow generated from the thermally produced particles after 
the decay of the curvaton flelds. In addition, if we consider the Affleck-Dine baryogenesis, 
Hb vanishes when the Affleck-Dine fleld starts to move after the beginning of the 0D epoch. 
If the baryon asymmetry is generated before the 0D epoch, on the contrary, fluctuation 
of the baryon density becomes negligibly small since baryogenesis occurs when the metric 
perturbations has not grown enough. Even if there is no primordial fluctuation in 

the baryonic sector, however, isocurvature fluctuation may arise in the baryonic component 



Fig.El 
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since density fluctuation of is inherited to tlie density fluctuation of tlie radiation. As 
we mentioned, Si becomes tlie entropy between components generated from 0j and those 
not from 0, and hence S'^^'^'*^ = —Y.if'y^sf'^^'' = |[\I/*-^'^'^^]rd2- In summary, the value of 
Kb depends on when the baryon asymmetry of the universe is generated: 

_ I : baryogenesis during/after 0D 
'^^ [9/2 : baryogenesis before 0D 

Next we consider the isocurvature fluctuation in the CDM sector. One important point 
is that, since the interactions of the moduli fields are very weak, it is natural to expect 
that the reheating temperature due to the decay of the moduli fields is very low. This fact 
has an important implication to the scenario where the lightest supersymmetric particle 
(LSP) becomes the CDM since the reheating temperature may be lower than the freeze- 
out temperature of the LSP. If so, they should be non-thermally produced to account for 
the CDM. In this case fluctuations in the curvaton fields may generate extra isocurvature 
fluctuations in the CDM sector as will be discussed below. Of course, there are other well- 
motivated candidates of the CDM, like the axion, with which the isocurvature fluctuation 
in the CDM sector vanishes.*^ 

Let us consider the most non-trivial case where the decay products of the moduli fields 
directly produce the LSP j2S|- In this case, the CDM density as well as the densities of 
the curvatons and radiation are determined by the following Boltzmann equations: 

dTi ~ — 

+ ?>Hn^^ = N^^iT^^n^^ - {v,^ia)n^n^^, (4.3) 

^"*^ + 3i7^^^ = -r^.n^., (4.4) 



dt 
dt 



+ ^Hp^^ = (m^^ - N^^im^)T^^n^^^ + m^{vreia)n^n^„ (4.5) 



where H = a/a, n^.^ is the number density of the LSP produced from the decay of (pi while 
= J2i total number density of the LSP, rriy. is the mass of the LSP, iV^^j is 

the averaged number of the LSP produced in the decay of one modulus particle 0j, and 
(vreicr) is the thermally averaged annihilation cross section of the LSP. (In our numerical 
calculation, we assume that the LSP is the neutral wino, as suggested by the simple 
anomaly mediated model. We expect that xx W^W~ is the dominant annihilation 
mode and we use 

271 (2 — xw) 

where xw = 'rn'wl'rn\ with mw being the W^-boson mass, and g2 is the gauge coupling 
constant of SU(2)l. Our conclusions are qualitatively unchanged even if the LSP is not 
wino.) 



axion starts its oscillation during/after </)D. 
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First, we consider the case with only one (real) modulus field. (For the moment, we 
drop the index i.) Number density of the relic LSP is sensitive to the parameters and 
(freicr). If A^^ -C 1, the pair annihilation of the LSP becomes ineffective and hence almost 
all the produced LSPs survive. Then, 

[n^U^r, - N^H^r, - (4-7) 
where tt,^ = p^/m^ being the number density of (p, leading to 



c 



2^3 

Parameterizing the decay rate of the modulus field as 



-iV.FfMy^^. (4.8) 



2-K M2 ' 



(4.9) 



we obtain 



^ ^ 5.8 X 10'^ GeV x iV,A, ( ) ( {^V^ (4.10) 

s * VlOO GeV/ VlOO TeV/ V10.75/ ^ ' 

Importantly, in this case, the resultant LSP density is insensitive to (freio")- On the 
contrary, if iV^ becomes large enough, the pair annihilation of the LSPs cannot be neglected 
and the number density of the LSP is given by 



H 



- r^v (4-11) 



Thus, the resultant LSP density is inversely proportional to (frcio"). 

Now we discuss the density fiuctuations in the multi-field case. As we discussed in 
the previous section, the metric perturbation is generated associated with the primordial 
fiuctuations of the moduli fields. In discussing the cosmic density fiuctuations, it is also 
important to understand the properties of the isocurvature fiuctuations. In the RD2 epoch, 
behaviors of the density and metric perturbations can be understood by using Eqs. (j3.19p 
and ()3.20|) . Isocurvature fiuctuation in the CDM sector is, on the other hand, obtained 
by solving Eqs. (lO) - 

In estimating the isocurvature fiuctuation in the CDM sector, it is important to note 
that 5*1 and 5*2 should be different to realize a non- vanishing value of Sc-y- If the ratio 
n^^/n^^ does not fluctuate, the system consisting of 0i and 02 can be regarded as a 
single fluid and hence the isocurvature fluctuation should vanish. This fact implies that, 
if — ^2 = 0, no isocurvature fluctuations can be induced as far as the components 
generated from the decay products of 0i and 02 are concerned. Thus, the isocurvature 
fluctuation should be proportional to 5*12 = — 5*2, and 
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Defining 



we obtain 



Notice that 



K 



(Ha) 



K 



(Ha) 



QiHA) 
\^i.HA) 



RD2 



'-'12 



dln[n^Jn^2 



(4.13) 



(4.14) 



(4.15) 



As one can see, Kc, which parameterizes the size of the correlated isocurvature fluc- 
tuation in the CDM sector, depends on various parameters. First, it depends on the 
parameters related to the properties of the moduli fields, in particular F^. and Nx,i- In 
addition, is sensitive to the initial amplitudes of the moduli fields on which depends. 
Furthermore, Kc depends on the transfer matrix T defined in Eq. (|3.15|) which parameter- 
izes how the primordial fluctuations in the moduli fields 6(f)A propagate to the fluctuations 
of the mass eigenstates 50^. In particular, T determines the ratio S^^'^'^^ s[2'^^\ Since the 
interactions of the moduli fields are suppressed by inverse powers of M*, the potential of 
is expected to be almost parabolic if the amplitude of </> is smaller than M*. Then, T 
becomes diagonal (in some basis) with a good approximation, hence 501 and 502 become 
uncorrelated. In this case, S / Su = ±|. 

In our study, we numerically solve Eqs. ()4.3p — (14. 5|) to obtain Kc given in Eq. ()4.14p . 
Although k[^'^'^^ has complicated dependence on various parameters, we can understand 
its behavior for several cases if T is given by Eq. ()3.16|) . Let us consider three extreme 
cases: 

1. F^^ < F^2, N^,i < 1 and [0i/02]inf^ > Max [m^^/m^^, N^^2m^2 / ^x^m^i 

The decay products of 02 is cosmologically negligible in this case, i.e., f-y^ ~ 1 and 
/ci ~ 1. Therefore, no isocurvature fluctuations are generated from the fluctuation 
of 01. On the contrary, 502 results in Sn^^/rij^^ ~ Sp^^/py^ ~ S'2. As a result. 



K 




The ratio of the metric perturbations is given as 



VI/W2) 



RD2 



fl2 
/71 



inf 



2/3 



2 r2/3 

inf m^^l ^2 



<1, 



(4.16) 



(4.17) 



where it should be noted that the masses of the moduli fields are of the same order 
of magnitude as mentioned in Sec. El Thus 502 does not contribute to Kcorr nor n^so 
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as long as N^^i ^ iV^,2, even though k^^'^^) jg finite (see Eq. ()3.30|) ). Thus we expect 
that the density fluctuation is purely adiabatic in this case. On the other hand, 
if [0i/02]inf ^^,2^^^/ ^xA^4'i is satisfied, Kjso can be large due to 502, so the 
isocurvature fluctuation becomes substantial. 



2. > r^2, N^j < 1 and [4>i/ 4>2]mf^ ^ N^,2m^J N^^iui^^: 



First we focus on the case of [0i/02]inf^ ^ y''^4>2^4>i/™ti^'l'2- ^^^^ case, most of 
the radiation and the CDM are generated from the decay products of 0i, and hence 
/-yj ~ 1 and ~ 1. Then we obtain the same result as in the previous case; 

^Wi)^0, ^ _9 ("^k _ (4.18) 



It is more interesting to consider the case of [0i/02]inf^~ ^m^2r^j/m^^r<^2. Then, 
most of the radiation comes from the decay of 02, while most of the CDM are 
generated from 0i, namely, f^^ — 1 ^ind ~ 1. Thus we expect large isocurvature 
fluctuation. 

3. = r<^2 and iV^,i < 1: 

In this case, (pi and 02 simultaneously decay. Assuming N^^i -C 1, almost all 
of the CDMs produced by the decay do not experience the pair annihilation and 
hence the resultant number of the LSP is approximately proportional to [iV-^^in^^ + 
^x,2'^02]H~r0- Using the fact that [m^jri,^^ + m^2'^02]-H'~r0 is fixed, we obtain 

2/,, (%i + %2a:)(l+y) ^ ^ 

where x = [n^^/n^-^]i^f and y = [m^^n^^/m^^n^^]i.[ii. k^^'^'^^ is obtained by replacing 
1^2 and x, y ^ 1/x, 1/y. Notice that the isocurvature fluctuation vanishes if 
Ny.^i/m^-^ = N-^fl/m^.^. This is from the fact that 0i and 02 can be regarded as a 
single fluid, since n^^/n^^ = ny.^/n-^^ when Ny.^i/m^-^ = N^^2/m^2- 

We numerically calculate K^orr and Hiso with the use of Eq. ()3.30|) . Contours of constant 
Kcorr and Kiso are shown in Figs. 121 and 0] Then, with Kcott and Hiso, the CMB angular 
power spectrum is calculated and we derive constraints on the curvaton model under 
consideration making use of the likelihood analysis as Fig. |21 In the numerical calculation 
and the discussion hereafter, we assume m^-^ = m^^ for simplicity. Also we take m-^ = 100 
GeV, and A^^ runs from 0.1 to 10. *^ 

We first consider the case with = 0. Even in this case, there are two scalar fields 0i 
and 02 and hence, in general, correlated mixture of the adiabatic and isocurvature fluctu- 
ations may arise. In order to discuss the effects of the correlated isocurvature fluctuation, 

"^^In our numerical calculation, we have tuned the mass of the moduli fields so that the pair annihilation 
cross section of the lightest neutralino realizes our canonical value of the relic density of the CDM. 
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10 %:., ' 




10 10 



Figure 3: Contours for the values of /tiso and Kcott in ri/r2 vs. [0i/02]inf plane for A^^_i = 
10~^, = 10~^ and Hb = 0. We also show the allowed area within 95 % CL. as the 
shaded (green) region. Here m^^ = m^^ is assumed. 



we calculate the likelihood variable and derive the constraint. In Figs. 01 and 01 we plot 
the contours for the probabilities obtained from the likelihood analysis. The results with 
A''^^! = N^^2 = 10~^ (see Fig. can be relatively easily understood. In this case, almost 
all the components are generated from 0i when ri/r2 ^ 1 and [0i/02]inf ^ 1 {i.e., the 
lower-right corner of the figures), while the decay products of 02 generates almost all the 
components when ri/r2 ^ 1 and [0i/02]inf 1 {i.e., the upper-left corner of the figures). 
In these cases, the situation is like the case with single scalar field and hence the CMB 
angular power spectrum almost agrees with the adiabatic one. Remarkably, in order to re- 
alize the "adiabatic-like" result, the ratios ri/r2 and [0i/02]inf do not have to be extremely 
large as can be seen in Fig. |H1 In the upper-right and lower-left corners of the figures, de- 
cay products of both 0i and 02 produce significant amount of the radiation and/or the 
CDM. Note that the second extreme case mentioned before cannot be seen in this figure, 
because A^^ i = = 10~^ are not small enough and hence the pair annihilation of the 
LSP is not negligible. We have found that l^corrl and /tiso are enhanced when the amount 
of the CDM from 0i is comparable to that from 02- In these cases, Scy becomes sizable 
and the shape of the CMB angular power spectrum becomes different from that from the 
adiabatic fluctuations. Indeed, the CMB angular power spectrum at lower multipoles is 
enhanced relative to that at higher multipoles due to this effect. Furthermore, as pointed 
out above, one can see that the isocurvature fluctuation vanishes along the horizontal axis 
with Fi = r2. Even if there exists large hierarchy between A^-j^_i and the constraint 
does not change much. In Fig. 01 we plot the result for the case N^^i ^ Ny.^2- 
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Figure 4: Same as Fig. El except A^^^i = 0.1 and = 10~^. 

When Kh 7^ 0, the constraint completely changes. As we mentioned, Kf, = 9/2 if 
the baryon asymmetry is generated before the epoch, assuming vanishing primordial 
fluctuation in the baryonic component. In this case, the correlated isocurvature fluctuation 
almost always becomes large, leading to unacceptably large Kcorr and Kiso irrespective of 
the choice of A^^^^i and Indeed, in this case, we have found that the predicted 

CMB angular power spectrum becomes consistent with the WMAP result in a very tiny 
parameter region. 

Thus we conclude that mild hierarchy in ri/r2 is necessary to realize the CMB angular 
power spectrum consistent with the WMAP results, in the case of vanishing baryonic 
isocurvature fluctuation. On the other hand, if there is nonzero baryonic isocurvature 
fluctuation, it is difficult to obtain the CMB angular power spectrum consistent with the 
observations. In order to have successful curvaton scenario, we need to suppress Kb as 
which would constrain the scenario of baryogenesis. 

5 Affleck-Dine Field as Curvaton 

Another class of possible curvaton fields in supersymmetric models includes F- and D-fiat 
directions. There are various fiat directions in the supersymmetric models, which are lifted 
only by the effect of the supersymmetry breaking. One important application of such a 
flat direction to cosmology is the Affleck-Dine baryogenesis where baryon asymmetry of 
the universe originates from the condensation of such a flat direction [24j. Thus, now we 
consider the case where the curvaton fleld is one of such flat directions. In particular, 
in this section, we study the possibility that the Affleck-Dine fleld plays the role of the 
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curvaton. In this case, the curvaton field is also responsible for the baryon asymmetry 
of the universe. Another case will be discussed in the next section. 

In the Affleck-Dine scenario, the baryon asymmetry of the universe is generated in the 
form of the coherent oscillation of the Affleck-Dine field. To make our discussion clearer, 
we adopt the following form of the potential 

V = m'M" + J^2i^'' + ^-^-)^ (5-1) 

where the scalar field has non-vanishing baryon number B^. Here, since is an F- 
and D-Hat direction, is induced by the effect of the supersymmetry breaking and is 
of order the gravitino mass.*^ In addition, A is a dimensionless coupling constant (which 
is assumed to be 0(1)) and n is an integer larger than 2. (We choose the coefficient of 
the higher- dimensional operator being real using the phase rotation of 0.) We expect that 
the baryon-number violating term is induced from a higher dimensional Kahler interaction 
and hence is proportional to the parameter m|.*^° As discussed in Sectional we assume 
that the potential is not affected by the Hubble-induced interactions. 
The equation of motion of the Affleck-Dine field is given by 

" dV 

+ 3i70+^ = O, (5.2) 



which leads to 



fdV \ 

fiB + SHub = iB^ — -(f) — h.c. , (5.3) 



where the baryon-number density in this case is given by ns = iB^{(f)*(j) — (p 0). Assuming 
non-vanishing initial amplitude, nonzero baryon number is generated when starts to 
oscillate. The baryon number density at that moment is estimated as 



dV 

BAm I ^ 




~^<^]^C.tSm(n^init), (5.4) 



where we have neglected coefficients of 0(1). Here, we have used the approximation that 
the effect of the baryon-number violating operator is so small that it can be treated as a 

"^^In this section we assume the gravity- mediation. If the gravitino mass TO3/2 is much larger than the 
soft mass m^, as in the case of anomaly- mediation ^JEIE], the right-hand sides of Eqs. (|5.4(l and (|5.f)|l 
are multiphed by a factor of im2/2/m^)^ . (In this case ^init ^ M^{m^/m^/2)'^^'""'~^^ is required to avoid a 
color /electromagnetic breaking minimum j25|.) However, the conclusion of this section does not change. 
We do not consider the case of gauge-mediation '5^, where in general a stable Q-ball is generated '27'. 
#iOrpj^jg is not the case if there is a non-rcnormalizable operator in the superpotential. In such a case, the 
higher order term in Eq. 15.1|l is proportional to m^, and Eqs. (|5.4I) and l|5.6() arc multiplied by (Af^/m^). 
However, the conclusion of this section docs not change, as long as the Affleck-Dine field plays the role of 
the curvaton. 
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perturbation. (Validity of such an approximation will be shown later.) As one can see, 
the baryon-number density becomes smaller as the initial amplitude of the Affleck-Dine 
field is suppressed. 

The energy density of the Affleck-Dine field is given by ~ m^l^p, where we neglected 
the effect of the higher dimensional operator assuming 0init is small enough. Thus, defining 



e = 



(5.5) 

(/)D 



the e parameter is given by 

^ ~ ^'t'^ ( ) sin(n6'init). (5.6) 




Importantly, when H ^ m^, hb and are both proportional to and hence the ratio 
TT's/Pcj) is a constant of time as far as F^. (In the case of Q-ball formation, and 
should be understood as the energy density and the decay rate of the Q-ball, respectively.) 

When the Affleck-Dine field oscillates fast enough with the parabolic potential, the 
motion of the Affleck-Dine field is described as 

= Ai(t) cos{m^t), (t)2 = A2(t) sm{m^t), (5.7) 

where we have used the phase rotation of the Affleck-Dine field to obtain the above form 
with Ai > A2. The amplitudes Ai{t) and A2{t) are proportional to a"^/^ when <^ 
H <^ m^. 

Evaluating the ratio nB/p4, at the time of the decay of the Affleck-Dine field (or the 
Q-ball), the resultant baryon-to-entropy ratio is estimated as 

— ~ e , (5.8) 

s 

where Trd2 is the reheating temperature due to the curvaton decay. In order not to spoil 
the success of the big-bang nucleosynthesis, Trd2 should be higher than about 1 MeV [2H] • 
In addition, is as large as the gravitino mass and is expected to be smaller than 
~ 1 TeV assuming the supersymmetry as a solution to the naturalness problem in the 
standard model. Consequently, to explain the presently observed value of the baryon- 
to-entropy ratio of ub/s ~ 0(10"^''), the e parameter should be much smaller than 1, 
that is, 10~^. This fact means that a large hierarchy between Ai and A2 is required: 
1^2/^1! ~ 10~^. Such a hierarchy is realized when, for example, the initial amplitude of 
the Affleck-Dine field is smaller than the suppression scale of the baryon-number violating 
higher dimensional operator (i.e., in our example, M^). The above fact tells us that the 
motion of the curvaton field is almost in the radial direction. In addition, the field 0^ 
and (f)0 (almost) correspond to the mass eigenstates 0i and 02? respectively, since the 
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scalar potential is well approximated by the parabolic one. As we will discuss, primordial 
fluctuations in 0^ and 00 give rise to different type of density fluctuations. 

Importantly, behaviors of the Affleck-Dine fleld strongly depends on the detailed shape 
of the potential. Even when the amplitude of is small enough so that the higher dimen- 
sional interaction in Eq. (jS.lj) can be disregarded, the scalar potential may deviate from 
the simple parabolic form once the renormalization group effects are taken into account: 

V = m|(M,) [l + inog(|0|VM2)] |0|2, (5.9) 

where we chose as a renormalization point. The parameter K is from the renormaliza- 
tion group effect. If the Yukawa interaction associated with is weak enough, the gauge 
interaction determines the scale dependence of and the parameter K then becomes 
negative. If the Yukawa interaction becomes strong, on the contrary, is more enhanced 
at higher energy scale and the parameter K becomes positive. Thus, if the Affleck-Dine 
mechanism works with the flat direction consisting of flrst and/or second generation MSSM 
particles, the parameter K is likely to be negative. If a third generation squark (in par- 
ticular, stop) or the up-type Higgs fleld is associated with the flat direction, however, K 
may become positive. 

Evolution of the coherent oscillation depends on the sign of K. If the potential of is 
flatter than the parabolic one, coherent oscillation of the scalar fleld may evolve into the 
non-topological solitonic objects called "Q-ball" (or in our case, "i?-baU") [22] • Properties 
of the density fluctuations differ for the cases with and without the Q-ball formation. 
Thus, we discuss these two cases separately. 



5.1 Case without Q-ball formation 

First, we consider the case where the Q-balls are not formed. This is the case when 
the parameter K is positive. In this case, it is expected that the reheating temperature 
after the decay of the Affleck-Dine fleld is relatively high, i.e., higher than the freeze-out 
temperature of the lightest neutralino which is assumed to be the CDM. In this case, 
the relic LSPs are thermally produced and relic density of the LSP (more precisely, nd s) 
becomes independent of the reheating temperature due to the decay of the Affleck-Dine 
fleld. Thus, in this case, no entropy fluctuation between the CDM and the radiation is 
generated. 

As we mentioned, when the amplitude of the Affleck-Dine fleld is much smaller than the 
suppression scale of the baryon-number violating higher dimensional interactions, the flelds 
(pr and 06) (almost) correspond to the mass eigenstates 0i and 02, respectively. Primordial 
fluctuations in 0^ and (pe result in different type of density fluctuations. Because of the 
large hierarchy between the amplitude of the flelds 0i ~ 0^ and 02 — 0e, almost all 
the photons in the RD2 epoch are generated from the decay products of 0i and hence 
/^^ ^ 1 while — 0- Furthermore, since the Affleck-Dine fleld mostly feels the parabolic 
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potential, the primordial entropy fluctuations are given by 



Si 



As a result, using Eq. ()3.20p . we obtain*^^ 



(5.10) 

J inf 



RD2 9 0;^.,' 



~0. (5.11) 

RD2 ^ ^ 



In addition, since S(j)r and 6(j)g provide fluctuations in 0init and 6, respectively, the entropy 
fluctuation between the baryon and radiation is also expected. Indeed, using the fact that 
the e parameter given in Eq. ()5.5|) is proportional to the resultant baryon-to-entropy ratio, 
entropy fluctuation in the baryonic sector is estimated: 

= (n - 2)^, ) =ncotin9,^,A (5.12) 

<Pinit rinit 



Since the fluctuations 6(j)r and 6(j)g are uncorrelated, the total CMB angular power 

-We) 
I ; 



spectrum is given by a linear combination of the two power spectra C/''^'-) and C/^"^"^ 
which are characterized by the following parameters: 



k(^^'-) ^ -^i^^-^, k(;1^«)^oo. (5.13) 

Thus, the CMB angular power spectrum associated with 6(j)0 is (almost) the same as the 
purely isocurvature result while C^^'^'''^ is from the correlated mixture of the adiabatic 
and isocurvature fluctuations. Relative size of C^'^'^'-) and Cl^^'^ depends on 9;^^^ bmce 
^m'^'^'' < 0' however, the acoustic peaks are always suppressed compared to the adiabatic 
case once one normalizes the Sachs- Wolfe tail. Since the angular power spectrum measured 
by the WMAP is well consistent with the adiabatic result, suppression of the acoustic peaks 
causes discrepancy between the theoretical prediction and the observations. To study this 
issue, in Fig. we plot the predicted CMB angular power spectrum for the case with n = A 
normalizing the Sachs- Wolfe tail. Here, we use the relation fl5.13p with cot ri^init = which 
makes the discrepancy smallest. As one can see, even in the case with cotn6'init = 0, the 
CMB angular power spectrum becomes extremely inconsistent with the WMAP results. 
Indeed, the total angular power spectrum is characterized by the following k parameters: 

9 (n - 2)fifc 9 VLi, 



4 ^\iyy^ 4 ^^7" 



TTT^V - 2)' + cot2 {ne,^^), (5.14) 



and, using n > 3 and adopting our canonical values of and fi^, we obtain Kcott ~ —0.39 
and l/tisol > 0.39, which is inconsistent with the observational results. (See Fig. El). 

'^-'^^In fact, associated with 54>0, the metric perturbation [^^'^'^''''] ^152 '^^y generated. However, it is 
suppressed by a factor of f^^'^^/ fjiSi ^ |A2/Ai|((5(/)e/(5(/)r) 10^'*, compared with [^^'''^'■^] r{i52' ^^^tric 
perturbation of this order does not significantly change the shape of the CMB power spectrum and we 
neglect such an effect. 
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Figure 5: The predicted CMB angular power spectra for the case where the Affleck-Dine 
field plays the role of the curvaton. Prom top to bottom, the WMAP data and the best-fit 
ACDM model, the case without the Q-ball formation, the case with the Q-ball formation 
[C = 1) which accidentally coincides with the previous case, the case with the Q-ball 
formation (C = 0), and the purely isocurvature case. We take n = 4, cotn6'init = and 
normalize the Sachs- Wolfe tail. 
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Therefore, the present case is excluded by the WMAP result. As we have shown, this 
is due to too large entropy fluctuation in the baryonic sector j^Oj. We should emphasize 
that the result obtained here is quite generic, independently of the detailed dynamics of 
the Affleck-Dine mechanism as far as the Affleck-Dine fleld is an F- and D-flat direction. 
Thus, we conclude that the Affleck-Dine field cannot play the role of the curvaton. In the 
next subsection, we will see that essentially the same conclusion is obtained for the case 
with the Q-ball formation. 



5.2 Case with Q-ball formation 

If the potential of the Affleck-Dine field is fiatter than the parabolic one, it is inevitable 
that the Q-balls are formed. This is due to the fact that, if < 0, instability band arises 
in the momentum space, which is given 

0</cJ,^, <2mJ|i^|, (5.15) 

where /Cphys is the physical momentum. In this case, once the Affleck-Dine field starts to 
oscillate, fiuctuation with the momentum within the instability band grows rapidly. 

Typical initial charge of the Q-ball generated in this process depends on the value 
of e. If e is close to 1, charge of the typical Q-ball is given by the total charge within 
the horizon at the time of the Q-ball formation. As a result, typical initial charge |Q|init 
is proportional to e. In this case, formation of the Q-ball with opposite sign of charge 
(so-called anti-Q-ball) is ineffective. If e becomes much smaller than 1, however, |Q|init 
becomes independent of e. In this case, almost same number of the Q-ball and anti-Q-ball 
are formed with the relation e = {uq — nQ)/{nQ + ng), where nq and hq are the number 
densities of the Q-ball and anti-Q-ball, respectively. Indeed, numerical lattice simulations 
have shown that 

|Q|,„,..4(!^yxj^ (5.16) 

\ I \ ec for e ^ ec 



where /3„ ~ 6 x 10-=^ and Cc ~ 0.01 ^*^^ 

In addition, the condition for the instability band (|5.15p provides the size of the Q-ball 

as 

~ 1^11/2 • (5-17) 



'^^^In fact, the numerical factor for the instability band depends on the value of e, but the order of 
magnitude does not change. 

'^^^In fact, the charges of produced Q-balls distribute with typical value given by Eq. H5.16|l . The precise 
estimation of the distribution is difficult since the Q-balls are produced through meta-stable state "/-ball" 
and it needs very long time simulations to find the distribution of the final stable Q-balls. However, 
the dynamics only depends on (/)init, *-e., which justifies the present analysis. 
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As long as the initial value of the Q-ball charge is larger than 0(10^^), evaporation of the 
Q-ball is ineffective and the universe is reheated by the decay of the Q-ball. The decay 
process of a single Q-ball is described by 

where A = 4:71 Rq is the surface area of the Q-ball (This relation is valid when 

the charge of the Q-ball is much larger than 1.) Thus, in our case, Q is (almost) a 
constant of time and hence the lifetime of the Q-ball with initial charge Qinit is given by 
Tq^ ~ IQinit/QI- Using Eq. ()5.16|) with e^ec, the decay temperature of the Q-ball, i.e., 
^RD2, is estimated as 



Trd2 ~ 1 MeV X 



li^l / \3/2 



0.03/ VlTeV/ 




One of important consequences of the Q-ball formation is that the RD2 epoch may be 
realized at much lower temperature compared to the case without Q-ball formation; for 
the Q-ball charge Qmit ~ 10^^ — 10^^, the RD2 epoch is realized at the temperature 
^RD2 ~ 1 MeV — 1 GeV, which is lower than the freeze-out temperature of the lightest 
neutralino (which is assumed to be the LSP). The other important consequence is that the 
decay rate of Q-balls, Tq, can be another source of the density fluctuations [HII. Note that 
Tq depends on (/)init, and that the decay process is different from the usual exponential 
decay. Hence, if 0init fluctuates spatially, Tq has spatial variation, leading to additional 
density fluctuations. 

Now we discuss the density fluctuation in the case with Q-ball formation. First, we 
consider the metric perturbation. Taking account of the fact that fluctuation of the de- 
cay rate generates adiabatic density perturbation [31], the metric perturbation has the 
following form: 

^WOl ^_(^ + c]^, \¥''^''>] ~0, (5.20) 

JrD2 V9 J 0i„it' L Jrd2 ' ^ ' 

where the second term in v]>('^9^'-) represents the contribution due to the varying decay 
rate, and the typical value of C is estimated to be order of unity. For our purpose, it 
is not necessary to know the precise value of C, but we crudely set C = 1. As long 
as C ~ 0(1), our conclusion does not change. Note that its sign is determined by the 
fact that the lifetime of Q-balls becomes longer for larger 0init, leading to more negative 
gravitational potential. Next we consider the isocurvature fluctuation in the baryonic 
and CDM sector. Since the density fluctuation due to the varying decay rate does not 
change entropy perturbations, we take Tq constant in the following. While the baryonic 
isocurvature fluctuation is calculated as in the previous subsection, the property of 6c can 
be different from the case without Q-ball formation. 
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To discuss the isocurvature fluctuation in the CDM sector, it is important to under- 
stand how the rehc density of the lightest neutralino {i.e., the CDM density) is determined. 
If Trd2 is lower than the freeze-out temperature of the LSP, the relic abundance of the LSP 
is non-thermally determined. The relic abundance depends on whether the annihilation 
rate of the LSP at the decay is larger or smaller than the expansion rate of the universe. 
If the annihilation rate is smaller than the expansion rate, (almost) all the LSPs produced 
by the decay of the Q-ball survive. In this case, no entropy fluctuation is expected be- 
tween the LSP {i.e., the CDM) and the photon produced by the Q-ball decay. If all the 
LSPs produced by the Q-ball survives, however, the relic density of the LSP is expected 
to become much larger than the critical density of the universe. This is because the Q-ball 
consists of the condensation of supersymmetric particles {i.e., the scalar quarks and/or 
scalar leptons) and hence the situation corresponds to the case discussed in Section HI with 
Ny. ~ 0.1 — 1. Indeed, from Eq. ()4.8|) . we can estimate the LSP density normalized by the 
entropy density of the universe: 

^ ~ 10-^GeV ( ^ ^ ^ ^ ^ """^ 



s 



0.1/ UMeV/ UOOGeV/ UTeV 



which is much larger than the observed value Pcrit/s — 3.6 x lO^^/i^GeV, with pcrit being 
the critical density of the universe. 

Thus, in the following we discuss the case where the annihilation rate of the LSP 
is larger than the expansion rate of the universe. In this case, the pair annihilation 
proceeds until the annihilation rate becomes comparable to the expansion rate and hence 
[nj^]H~rQ ~ Tq/ (^^reio")- In the RD2 epoch or later, ratio of to is conserved and hence 
we obtain 

^ ~ — r, 5.22 

where we have used the relation Tq ~ T^j)2/^*- We assume that Tq and (froicr) satisfy 
the relation such that the above ratio becomes consistent with the currently observed dark 
matter density ^3] . 

As shown in Eq. ()5.19|) . Trd2 depends on 0init and hence the ratio ny./n^ fluctuates if the 
Affleck-Dine field has primordial fluctuation. Using Trd2 oc ^J^it, we obtain n^jn^ oc 0init- 
Then, 

5W^) = i^, ^We) = o, (5.23) 

0init 

where we have used the relation S'c^ = 5{nc/n^) / {ric/n-^). In this case, the fluctuation in 
the radial direction, 50^, also generates the isocurvature fluctuation in the CDM sector. 
Notice that the fluctuation in the phase direction 5(f)e does not affect the isocurvature 
fluctuation in the CDM since e -C 1. 
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In the case of the Q-ball formation, significant amount of the isocurvature fluctuations 
are generated associated with dipr and 5(j)0. Adding the isocurvature fluctuations in the 
baryon and the CDM, we obtain 



Importantly, n^^'''' is again negative and k!:^'^^ ^ 1. Thus, 50r results in the density fluc- 
tuations with correlated mixture of the adiabatic and isocurvature fluctuations while 5(f)e 
becomes (almost) purely isocurvature fluctuations. The k parameters, which characterize 
the total CMB power spectrum, are thus given by 



For the total CMB angular power spectrum, if one normalizes the Sachs- Wolfe tail, 
acoustic peaks are extremely suppressed compared to the purely adiabatic case. In Fig. El 
the CMB angular power spectrum with the relation Eq. ()5.25j) is also plotted in the case 
of cotn^init = and C = and 1. In this case, again, the shape of the total CMB angular 
power spectrum becomes inconsistent with the observations. 

To summarize, the Affleck-Dine field cannot be the curvaton in both of the cases with 
and without Q-ball formation because the entropy density fluctuation in baryonic and/or 
CDM sector is too large. In other words, if Affleck-Dine fleld (or Q-ball) dominates the 
energy density of the universe, its primordial fluctuation must be suppressed, and the 
dominant density fluctuation must be generated by another source, like inflaton. 

6 F- and D-Flat Direction as Curvaton 

In this section, we consider the case where an F- and D-flat direction plays the role of 
curvaton which is not related to the mechanism of baryogenesis. Even in this case, Q- 
ball may or may not be produced depending on the shape of the curvaton potential; in 
particular, if the curvaton potential is flatter than parabolic (which may be due to the 
renormalization effects), Q-ball can be produced once the curvaton starts to oscillate. 
As can be expected from the discussion given in the previous section, properties of the 
cosmological density fluctuations strongly depend on how the curvaton fleld evolves. 

First, let us consider the simplest possibility, namely, the case without the Q-ball 
formation. In this case, Trd2 is expected to be higher than the freeze out temperature of 
the lightest supersymmetric particle and we assume that the CDM is thermally produced 

'^^'^In fact, the primordial fluctuation of the Affleck-Dine field is suppressed if there are Hubble-induced 
terms in the potential and/or the initial amplitude is large enough. 




~ OO. 



(5.24) 




(5.25) 
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lightest neutralino. If Trd2 is high enough, baryogenesis may be possible with the decay 
products of 0. For example, the Fukugita-Yanagida mechanism jSH] using the decay of 
thermally produced right-handed neutrinos may work if Trd2 ~ 10^"^° GeV p]. Another 
possibility may be electroweak baryogenesis ^ if Trd2 ^ 100 GeV. 

In this simplest case, no isocurvature fluctuation between the non-relativistic compo- 
nents and the radiation is generated by 50r nor dcpe {i-e., Km = 0). Thus, in the RD2 
epoch, the only source of the cosmic density fluctuations is the metric perturbation given 
in Eq. ()3.20|) . Since no isocurvature fluctuation is generated, Ci^"^''^ and cf'^"^ are both 
from purely adiabatic density fluctuations and hence cf'^'''' oc Ci^'^''\ Thus, even though 
\^{s<t>r) g^j-^^ \^iS(l)e) g^Ye free parameters, it does not affect the shape of the total CMB angu- 
lar power spectrum and Ci is the same as the result with purely-adiabatic scale-invariant 
primordial density fluctuations. (Concrete models for the flat-direction curvaton is found 
inRef. EH].) 

As discussed in the previous section, if the scalar potential of the flat direction field 
is flatter than the parabolic one, the Q-balls are formed. If the initial charge of the Q- 
ball is small enough, it decays (or evaporates) with a temperature above the freeze-out 
temperature of the LSP. This case is reduced to the one discussed just above, where 
the purely-adiabatic density fluctuation is obtained. If the charge of the Q-ball is large, 
however, it decays after the thermal pair-annihilation of the LSPs are frozen out. Thus, 
one can see that the entropy fluctuation is generated in the CDM sector: 



Si?''^ = T^, S^y = 0, (6.1] 

-^init ' 



which leads to 



+ ^^'^' = 0, (6.2) 

where C ~ 0(1) represents the contribution due to the varying decay constant. Notice that 



there is no baryonic entropy 
asymmetry.*-*^^ In this case. 



iuctuation, since the curvaton has nothing to do with baryon 



and hence 6(j)0 does not play any significant 



role in generating the cosmic density fluctuations. Thus, with the best-fit values of il^ and 
Qm obtained from the WMAP experiment, /tcorr — —1-9, —0.57 for C = 0, 1 respectively, 
which are clearly excluded by the observation. (See Fig.|21). Therefore, again, the curvaton 
with Q-ball formation which decays after the freeze-out of the LSP is rejected. 



7 Right-Handed Sneutrino as Curvaton 

Recent neutrino-oscillation experiments strongly suggest tiny but non-vanishing masses 
of the neutrinos. Such small masses can be naturally explained by introducing the right- 

'^-'^^However, if baryogenesis takes place before the Q-ball dominated universe, large isocurvature fluctua- 
tion is induced in the baryon sector, i.e., Kb = |. Here we assume that the baryon asymmetry is somehow 
generated after the beginning of Q-ball-dominated universe. 
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handed neutrinos P^ . 

If the right-handed neutrino exists in supersymmetric models, its superpartner may 
have non-vanishing primordial amplitude. In particular, such a case has been studied 
related to the possibility of generating the lepton number (and baryon number) asymmetry 
of the universe from the decay product of the right-handed sneutrino [10] • In this case, 
the right-handed sneutrino can be an another candidate for the curvaton |1T]. 

The scenario with the right-handed sneutrino condensation is almost the same as the 
previous cases except for the origin of the baryon asymmetry of the universe. During 
inflation, the right-handed sneutrino is assumed to have a non-vanishing amplitude. After 
inflation, the right-handed sneutrino starts to oscillate when the expansion rate of the 
universe becomes comparable to the sneutrino mass m^. If the decay of the oscillating 
right-handed sneutrino takes place much later than the infiaton decay, the sneutrino dom- 
inates the universe. Then, the sneutrino decays when the expansion rate of the universe 
becomes comparable to its decay rate. At the time of the decay, lepton number asymmetry 
can be generated provided that CP violation exists in the neutrino sector. Once the lepton 
number asymmetry is generated, it can be converted to the baryon number asymmetry 
due to the sphaleron process P?] . 

In this scenario, quantum fluctuations that the sneutrino acquires during inflation 
accounts for the adiabatic density fluctuations and hence the right-handed sneutrino plays 
the role of the curvaton. Here, baryon asymmetry of the universe is generated from the 
decay products of the right-handed sneutrino and hence no isocurvature fluctuation is 
expected in the baryonic sector. In addition, the decay rate of the right-handed sneutrino 
is high enough so that the LSPs (and other superpartners of the standard-model particles) 
are thermally produced. Then there is no isocurvature fluctuation in the CDM sector if 
the LSP becomes the CDM. 

If this is the case, the primordial density fluctuations become purely adiabatic and 
hence the resultant CMB angular power spectrum becomes well consistent with the results 
of the WMAP experiment. Thus, we conclude that the right-handed sneutrino is one of 
the good and well-motivated candidates for the curvaton. 

8 Conclusions and Discussion 

We have studied several candidates for the curvaton in the supersymmetric framework. 
Since scalar fields appearing in supersymmetric models are inevitably complex, it is nec- 
essary to consider cases with multi-curvaton fields. In this case, primordial fluctuations 
in the curvatons induce isocurvature (entropy) fluctuations as well as the adiabatic ones. 
Such isocurvature fluctuations may affect the CMB angular power spectrum and the recent 
WMAP experiment sets stringent constraints on those models. 

One potential problem of the curvaton scenario in the supersymmetric models is to 
suppress the baryonic isocurvature fluctuations. As shown in the present paper, isocur- 
vature fluctuation in the baryonic sector may arise in various situations. In many cases. 
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such a baryonic isocurvature becomes too large to be consistent with the WMAP result. 

Besides baryonic isocurvature fluctuations, we have found that the CDM isocurvature 
fluctuations are important for moduli or Q-ball curvaton where the abundance of the 
dark matter [i.e., the LSP) is determined non-thermally and hence its fluctuations are 
isocurvature with correlation to the curvaton fluctuations. 

If the cosmological moduli fields play the role of the curvaton, it is possible that the 
reheating temperature due to the decay of the moduli fields is lower than the freeze-out 
temperature of the lightest neutralino. In this case, the CDM is non-thermally produced if 
the lightest neutralino is the CDM, and the entropy fluctuation between the two curvatons 
{i.e., the real and imaginary parts of the curvaton) in general becomes the isocurvature 
fluctuation in the CDM sector. This fact excludes some part of the parameter space but, 
as we have seen, slight hierarchies between the initial amplitudes and/or the decay rates of 
two curvaton fields are enough to make the resultant CMB anisotropy consistent with the 
WMAP result. Of course, if the lightest neutralino is not the CDM, then we may evade 
the constraint from the isocurvature fluctuation in the CDM sector. An example for such 
CDM candidates is axion. 

If the Affleck-Dine field becomes the curvaton, on the contrary, too large baryonic 
isocurvature fluctuation is inevitably induced. As we have seen, the acoustic peaks are 
always too much suppressed to be consistent with the observations if the Affleck-Dine field 
plays the role of the curvaton. Thus, the Affleck-Dine field cannot play the role of the 
curvaton. It is also notable that, in the case with Q-ball formation, extra isocurvature 
fluctuation may arise if the reheating temperature due to the Q-ball decay is lower than 
the freeze-out temperature of the lightest neutralino (which is assumed to be the CDM). 

For the case where the F- and D-flat direction (without baryon number) becomes 
the curvaton, isocurvature fluctuation in the CDM sector may also arise if the Q-ball is 
formed. In this case, the resultant CMB angular power spectrum is too much affected by 
the isocurvature fluctuation. Without the Q-ball formation, however, the F- and D-flat 
direction (but not the Affleck-Dine fleld) can be a candidate of the curvaton. 

Another good and well-motivated candidate of the curvaton is the right-handed sneu- 
trino. If the right-handed sneutrino acquires the primordial quantum fluctuations during 
inflation, it becomes purely adiabatic density fluctuations and the resultant CMB angular 
power spectrum can be well consistent with the WMAP result. 

In summary, for a viable scenario of the curvaton, the isocurvature contribution should 
be so small that the CMB angular power spectrum becomes consistent with the WMAP 
result. In the supersymmetric case, this fact provides stringent constraints on the curvaton 
scenario. Thus, it is necessary to look for curvaton candidates which do not generate the 
isocurvature fluctuations; some of them are already found in the present study. 
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